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ABSTRACT 


For  seven  earthquakes  a  comparison  is  made  between  the 
TFO  long-period  horizontal  array  and  the  sum  of  the  horizontal 
strain  and  pendulum  instruments  at  QCAZ.  The  performance  of 
both  systems  was  severely  limited  by  noise  of  non-seismic 
origin. *Within  these  limits,  the  performance  of  the  two  systems 
is  about  equal.  This  report  is  preliminary  in  the  sense  that 
the  installation  of  instruments  at  QCAZ  was  not  complete  at  the 
time  the  data  were  obtained.  A  lower  level  of  non-seismic  noise 
is  to  be  expected  when  installation  is  complete. 

Detailed  analysis  of  one  Chiapas,  Mexico,  earthquake  demon¬ 
strates  that  strain-pendulum  combinations  can  be  used  to 
enhance  and  identify  both  body  and  surface  waves. 


*Detail  comparison  is  made  for  one  earthquake  between  QCAZ 
and  the  beamed  vertical  and  horizontal  array  at  TFO. 
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INTRODUCTION  AND  THEORY 


The  new  Queen  Creek,  Arizona,  (QCA2)  strain  seismograph 
system  is  about  125  km  south  of  the  Tonto  Forest  Observatory 
(TFO) .  Figure  1  is  an  index  map  of  Arizona  which  shows  the 
relation  of  QCAZ  to  the  TFO  long  period  array.  The  strain 
seismograph  installation  at  QCAZ  is  described  by  Fix  and 
Sherwin  (1970). 

Seven  earthquakes  were  selected  for  analysis  whose  arrival 
azimuths  at  QCAZ  would  be  near  the  installed  azimuths  of  the 
horizontal  strain  seismographs,  55° (235°)  and  325° (145°). 

Figure  2  is  an  equal  azimuth  world  map  centered  at  QCAZ  on 
which  are  plotted  the  locations  of  the  seven  earthquakes  used 
in  this  study.  USCSGS  data  for  these  events  are  listed  in 
Appendix  I. 

Benioff  (1962)  described  the  use  of  combinations  of  hori¬ 
zontal  strain  and  pendulum  seismometers  to  enhance  body  and 
surface  waves.  The  response  of  a  horizontal  strain  seismometer 
to  longitudinal  plane  waves  contains  a  factor  cos2a,  where  a 
is  the  angle  between  the  strain  rod  and  the  direction  of 
propagation  of  the  wave.  The  response  of  a  horizontal  pendulum 
contains  the  factor  cos  a  (Benioff,  1935).  Figure  3  shows  the  azi¬ 
muthal  response  of  a  horizontal  strain  seismometer  and  a  hori¬ 
zontal  pendulum  to  longitudinal  waves  (or  to  the  horizontal 
component  of  Rayleigh  waves).  When  the  outputs  are  summed,  the 
combined  response  is  as  shown  in  Figure  4.  To  interpret  this 
plot  it  is  important  to  remember  that  the  azimuth  to  an  event 
is  180°  opposite  to  its  direction  of  propagation.  In  our 
applications  both  the  sum  and  difference  of  the  outputs  have 
usually  been  calculated.  For  the  55°  (northeast)  azimuth  pair 
the  sum  is  a  beam  to  the  southwest  and  the  difference  (strain 
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minus  pendulum)  is  a  beam  to  the  northeast.  In  order  to  sum  the 
outputs,  it  is  necessary  that  response  of  the  instruments  be 
matched.  This  is  commonly  accomplished  in  the  recording  process. 

The  response  of  the  strain  seismograph  to  transverse  waves 
(SH  or  Love  waves)  is  proportional  to  (sin  os)  •  (cos  a),  where  a 
is  the  angle  between  the  direction  of  propagation  of  the  earth¬ 
quake  and  the  installed  azimuth  of  the  strain  rod  (Benioff,  1935). 
Figure  5  shows  the  response  of  a  horizontal  strain  seismograph 
to  transverse  waves.  An  orthogonal  pair  of  strain  seismographs 
can  be  used  to  enhance  transverse  waves,  providing  the  direction 
of  propagation  is  between  about  30°  to  60°  from  the  installed 
azimuth:  this  is  accomplished  by  computing  the  sum  (to  cancel) 
and  difference  (to  enhance)  of  the  wave. 

Transverse  phases  whose  direction  of  propagation  is  within 
20°  of  the  installed  azimuth  can  be  recognized  by  the  fact  that 
they  are  recorded  by  the  transverse  pendulum  but  not  by  either 
of  the  orthogonal  strain  instruments.  Strain  seismographs  can 
not  be  rotated  to  produce  radial  and  transverse  components  as 
are  pendulum  seismographs.  The  reason  for  this  is  that  the 
transformation  equation  to  a  new  coordinate  system  contains  a 
shear  component  of  strain  which  is  not  recorded  by  an  orthogonal 
pair  of  instruments.  A  third  strain  seismograph  at,  e.g. ,  45° 
to  the  orthogonal  pair  can  be  used  to  record  the  shear  component, 
providing  that  the  usual  plane  strain  assumptions  are  applicable 
(Jaeger,  1956). 

In  comparing  the  TFO  long-period  horizontal  array  to  the 
QCAZ  strain  installation  it  was  necessary  to  eliminate  some 
elements  from  the  TFO  array.  One  instrument,  at  site  2,  was 
inoperative  during  the  period  19  to  25  May  1970,  when  the 
comparison  was  made,  so  this  site  has  thus  been  eliminated  from 
the  computations.  There  also  exists  non-stationary,  impulsive, 
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low-frequency  (T  -  60  secs)  noise  on  the  horizontal  long-period 
instruments,  which  was  large  enough  to  contaminate  the  comparison. 
Channels  on  which  this  noise  occurred  have  in  some  cases  been 
eliminated.  A  similar  type  of  noise  occurred  on  the  325°  azimuth 
strain  instrument  at  QCAZ.  Commonly  the  effect  of  this  noise 
was  eliminated  by  filtering  with  a  15  to  50  second  pass-band  and 
by  selecting  samples  which  as  far  as  possible  did  not  include 
the  noise.  Examples  of  the  noise  described  are  marked  with  a 
dashed  line  block  on  several  figures  in  the  results  section  of 
this  report. 

Figures  6  and  7  are  a  comparison  of  the  theoretical  hori¬ 
zontal  beam  response  pattern  for  TFO  and  QCAZ.  A  235°  azimuth 
beam  has  been  chosen,  corresponding  to  the  55°  azimuth  strain 
installation  at  QCAZ.  This  beam  is  toward  the  South  Pacific 
region  from  which  a  number  of  earthquakes  are  recorded.  In  both 
cases  for  TFO  a  20  second  period,  3.0  km/sec  wave  has  been 
chosen  to  determine  the  response.  At  QCAZ,  it  is  assumed  that 
the  local  phase  velocity  is  constant  with  period,  and  that  the 
strain  and  pendulum  instruments  have  been  matched  in  amplitude. 
Figure  6  is  the  comparative  beam  response  with  one  element 
(No.  2)  of  the  TFO  array  eliminated. 

Figure  7  is  the  beam  response  with  three  elements  (Nos.  2, 

5,  and  6)  eliminated.  The  response  of  the  TFO  array  includes 
the  cos  ol  response  of  the  horizontal  instruments  steered  for  a 
235°  azimuth  event.  The  comparative  3  dB  apertures  of  the  beams 
are: 

QCAZ  64° 

TFO  (six  elements)  48° 

TFO  (four  elements)  51° 

At  TFO,  the  long  period  instruments  are  used  in  this 
analysis.  Radial  and  transverse  components  are  computed 
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from  the  horizontal  instruments  for  a  given  earthquake 
azimuth.  Figure  8  is  an  example  of  the  rotation  output  for  three 
sites.  The  radial  and  transverse  seismograms  are  designated  R 
and  T  respectively.  In  this  case  the  Love  wave  is  well-defined 
on  the  transverse  instrument  at  about  01  10  and  the  Rayleigh 
wave  about  three  minutes  later  on  the  radial  instrument.  This 
visual  test  was  used  to  verify  that  the  proper  relative  weights 
(calibrations)  were  used  when  the  instruments  were  rotated.  The 
phased  sum  (PS)  and  the  unphased  sum  (US)  were  computed  for  the 
available  elements  of  the  TFO  array.  A  surface  wave  phase  velocity 
of  3.5  km/sec  was  used  to  compute  the  phased  sum.  Because  of 
the  lower  phase  velocity,  body  waves  are  enhanced  on  the 
unphased  sum  and  Rayleigh  waves  on  the  phased  sum. 

At  QCAZ  the  weighted  sum  and  difference  of  a  horizontal 
pendulum  and  strain  seismograph  have  been  computed  for  comparison 
with  the  TFO  phased  sum.  The  weighting  factor  includes  local 
phase  velocity  and  calibration,  and  has  been  adjusted  to  obtain 
maximum  cancellation  of  the  back  beam.  Other  procedures  used  on 
the  QCAZ  data  are  described  in  the  results  section. 

The  data  from  QCAZ  were  obtained  during  the  latter  part  of 
May  1970,  prior  to  the  installation  of  the  vertical  strain 
instruments.  Although  one  can  in  theory  obtain  an  equivalent 
vertical  strain  by  adding  the  outputs  from  two  orthogonal 
horizontal  instruments,  this  procedure  has  not  been  used  in 
this  report.  There  was  considerable  installation  activity  in 
the  mine  at  the  time,  which  apparently  reduced  data  quality. 

The  above  discussion  of  beam  comparison,  and  QCAZ  vs  TFO 
has  its  first  application  to  surface  waves.  The  combination  of 
strain  and  pendulum  instruments  can  also  be  used  to  recognize 
and  enhance  body  waves.  The  formulae  for  a  homogeneous  half-space 
relating  strain  and  pendulum  displacement  involve  the  angle  of 
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incidence  of  the  particular  body  wave,  the  elastic  parameters  of 
the  medium,  the  horizontal  component  of  phase  velocity,  and  the 
wavenumber  (Appendix  II).  Experience  obtained  in  the  application 
of  these  formulae  at  various  strain  sites  is  that  they  are 
principally  useful  to  establish  phase  relations. 

Romney  (1964)  derives  the  relation 

9u  3u  A  +  2y  3u 

x  y  «- 

3x  3y  A  32 


which  says  that  the  sum  of  the  horizontal  strains  is  equal  to  a 
constant  times  the  vertical  strain.  This  equation,  which  also 
assumes  homogeneity,  can  be  used  to  estimate  the  ratio  (A  +  2p)/A. 
Applying  this  equation  to  WMO  data,  we  found  that  A  =  2y 
is  a  better  approximation  than  the  usual  assumption  A  *  y.  At 
the  present  time  such  data  are  not  available  at  QCAZ. 

Strain-pendulum  combinations  can  also  be  used  to  determine 
local  phase  velocity  for  a  dispersed  surface  wave  train.  The 
technique  consists  of  measuring  the  relative  amplitude  and  period 
of  each  half-cycle.  Amplitudes  are  measured  only  at  the  peaks  and 
troughs  of  the  waves  to  avoid  the  inaccuracy  involved  in  comparing 
the  ratios  of  small  numbers.  The  mathematical  justification  for 
this  procedure  is  included  in  the  Discussion  of  Results.  Mikumo 
and  Aki  (1964)  determined  local  phase  velocity  by  this  method, 
using  the  strain  installation  at  Isabella,  California. 
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DISCUSSION  OF  RESULTS 


QCAZ-TFO  comparison 

Seven  earthquakes  recorded  at  QCAZ  and  TFO  have  been  selected 
to  compare  beam  response.  At  QCAZ  the  beam  is  the  sum  differ¬ 
ence  of  the  horizontal  long-period  strain  and  pendulum  seismo¬ 
graphs.  At  TFO  the  beam  is  the  phased  sum  of  the  steered  long- 
period  seismographs.  USC$GS  data  for  the  earthquakes  are  listed 
in  Appendix  I. 

The  ratio  of  signal  amplitude  to  rms  amplitude  of  the  noise 
is  computed  for  one  of  the  comparative  beam  response  studies;  the 
high  level  of  non-stationary  noise  gave  unsatisfactory  results 
for  other  examples  tried. 

Figure  9  compares  the  QCAZ  and  TFO  beams  for  an  earthquake 
off  the  coast  of  Chiapas,  Mexico, on  19  May  1970.  P,  S  and  the 
Rayleigh  wave  (LR)  have  been  marked.  For  this  earthquake, 
southeast  of  QCAZ,  the  sum  of  the  strain  and  pendulum 
seismographs  should  enhance  these  phases.  Channels  1,  3 
and  4  at  TFO  are  the  same  data  used  in  the  example  describing 
the  rotation  procedure  (Figure  8).  P  and  S  are  marked  on  the 
unphased  sum  at  TFO.  At  the  scale  of  these  plots,  the  travel 
time  differences  across  TFO  for  P  and  S  are  negligible.  Examples 
of  the  non-stationary  noise  are  enclosed  in  a  dashed-line  block. 

Figure  10  shows  a  19  May  1970  earthquake  from  the  Fiji  Islands 
region;  only  the  Rayleigh  wave  has  been  identified.  The  path  for 
this  earthquake  is  largely  oceanic.  The  Chiapas  earthquake  ( 

(Figure  9)  has  a  continental  path.  The  surface  waves  from  a 
Greenland  Sea  earthquake  occur  prior  to  the  Fiji  earthquake 
surface  waves.  The  azimuth  of  this  event  falls  outside  the  beam 
at  QCAZ;  thus  the  sum  and  difference  for  this  segment  of  data 
has  not  been  computed. 
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Figures  11  and  12  consist  of  two  presentations  of  the  QCAZ 
and  TFO  data  for  overlapping  signals  from  three  earthquakes. 

This  is  the  type  of  data  one  wishes  to  have  in  order  to  further 
evaluate  the  two  systems;  the  search  for  better  examples  is 
continuing.  There  are  two  Chiapas  earthquakes;  a  small  one 
followed  about  six  minutes  later  by  the  magnitude  4.6  earthquake 
listed  in  Appendix  I.  The  third  earthquake  is  from  the  New 
Hebrides  region.  Travel  time  table  analysis  shows  that  several 
body  wave  phases  from  the  New  Hebrides  earthquake  should  occur 
at  nearly  the  same  time  as  the  surface  waves  from  the  Chiapas 
earthquakes.  In  Figure  11,  both  QCAZ  and  TFO  are  beamed  to  the 
southeast.  The  only  well-defined  phase  is  the  Rayleigh  wave 
from  the  second  (larger)  Chiapas  earthquake;  other  tentative 
identifications  have  been  made. 

Figure  12  shows  the  same  data  as  Figure  11  except  that  here 
the  55°  azimuth  strain  and  pendulum  at  QCAZ  have  been  used.  TFO 
has  again  been  beamed  to  the  southeast,  but  the  horizontal  instru¬ 
ments  have  been  rotated  to  obtain  a  radial  component  directed  to 
the  southwest  (New  Hebrides).  At  QCAZ  the  Love  waves  from  Mexico 
should  be  well  recorded  only  on  the  55°  pendulum  instrument* 

Both  the  strain  and  pendulum  instruments  should  record 
body  wave  phases  from  the  New  Hebrides  earthquake,  and 
these  should  be  enhanced  on  the  sum  channel.  Phases  marked 
on  the  QCAZ  sum  channel  (S55  plus  P55)  occur  at  the  time  computed 
from  the  Jeff reys-Bullen  travel- time  curves.  The  TFO  sample 
includes  body  wave  phases  from  the  New  Hebrides  and  the  Love 
waves  from  Chiapas.  Body  wave  phases  should  be  enhanced  on  the 
unphased  sum.  In  computing  the  phased  sum  travel  time  delays 
for  surface  waves  from  the  southeast  have  been  used.  Thus  the 
Love  waves  from  Chiapas  should  be  enhanced  on  the  phased  sum, 
because  the  transverse  Love  wave  component  is  approximately 
aligned  with  the  radial  instruments  steered  southwest. 
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Figure  13  is  the  QCAZ-TFO  comparison  for  a  Kermadec  Islands 
earthquake  on  25  May  1970.  A  number  of  body  wave  phases  have 
been  identified  in  addition  to  the  Rayleigh  wave. 

As  a  means  of  comparing  the  opposite  azimuth  cancellation 
obtained  for  the  QCAZ  data,  the  earthquake  data  of  Figure  13 
(Kermadec  Islands)  has  been  steered  toward  the  opposite  azimuth. 

The  phased  sum  of  Figure  14  should  be  compared  with  the  difference 
channel  (S55  -  P55)  of  Figure  13. 

Figure  15  is  the  QCAZ-TFO  comparison  for  an  earthquake  south 
of  Panama.  In  this  case  there  is  a  substantial  character  difference 
in  the  Rayleigh  wave  between  QCAZ  and  TFO. 

In  order  to  obtain  a  numerical  comparison  of  the  effective¬ 
ness  of  the  beam-forming  at  QCAZ  and  TFO,  a  signal  and  a  noise 
sample  have  been  selected.  The  signal  data  used  are  from  the 
Fiji  earthquake  of  19  May  1970  (Figure  10).  The  non-stationary 
noise  on  the  strain  instrument  at  QCAZ  and  on  the  horizontal 
pendulums  at  TFO  restricted  the  choice  of  noise  samples.  The 
noise  sample  selected  consists  of  500  seconds,  beginning 
1050  seconds  after  the  beginning  of  the  data  sample  of  Figure 
10.  It  thus  occurs  after  the  Rayleigh  wave  from  the  Greenland 
Sea  earthquake  and  prior  to  the  noise  on  the  TFO  5R  instrument. 
Figures  16  and  17  are  plots  of  the  noise  samples  used  as 
recorded  on  the  horizontal  instruments.  Figure  16  is  selected 
from  TFO  and  Figure  17  is  the  corresponding  sample  selected 
at  QCAZ.  At  QCAZ  the  gains  have  been  adjusted  so  that  the 
Rayleigh  wave  amplitudes  are  equal  on  the  strain  and  pendulum 
seismographs.  This  amounts  to  compensating  for  the  average 
local  phase  velocity  (Figure  10).  Visual  inspection  of  the 
noise  samples  indicates  that  the  sample  selection  has  fulfilled 
its  purpose  by  avoiding  the  Rayleigh  wave  from  the  Greenland  Sea 
earthquake  and  the  non-stationary  noise  which  occurs  at  both  sites. 
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The  sample  has  been  filtered  with  a  15  to  50  second  band-pass 
filter. 

The  following  signal  to  rms  noise  ratios  for  the  Fiji 
earthquake  were  obtained  for  sum  and  difference  of  the  hori¬ 
zontal  instruments  at  QCAZ  and  the  horizontal  and  vertical 
instruments  at  TFO.  Figure  18  is  a  plot  of  the  beamed  vertical 
long  period  instruments. 


QCAZ  Data 

S/rms  db 

S55L  (Horizontal 

s  train 

55°  azimuth) 

10. 

16  20.9 

P55L  (Horizontal  pendulum  55°  azimuth) 

11. 

0  20.9 

SUM  (S55L  +  P55L) 

15. 

0  23.5 

DIF  (S55L-P55L) 

4. 

0  12.0 

TFO  Data 

Radial 

Vertical 

Channel  No. 

S/rms 

db 

S/rms 

db 

1 

10.35 

20.3 

8.64 

18.7 

3 

8.19 

18.3 

12.28 

21.8 

4 

5.87 

15.4 

10.49 

20.7 

5 

4.66 

13.4 

4.57 

13.2 

6 

-- 

-- 

5.17 

14.3 

7 

4.49 

13.0 

8.41 

18.5 

Mean 

6.71 

16.5 

8.26 

18.3 

Unphased  Sum 

3.14 

9.9 

4.40 

12.9 

Phased  Sum 

6.45 

16.2 

17.71 

25.0 

Thus  the  SUM  channel  at  QCAZ  has  1.5  db  less  signal  to  noise 
improvement  than  the  phased  sum  of  the  vertical  instruments  at 
TFO  and  7.3  db  greater  signal  to  noise  improvement  than  the 
phased  sum  of  the  radial  instruments.  At  QCAZ  the  SUM  channel 
is  2.6  db  better  than  the  radial  pendulum  (S55L).  It  seems 
likely  that  the  results  for  the  TFO  radial  instruments  is 
anomalous.  The  S/rms  ratio  of  the  phased  sum  is  0.3  db  less 
than  the  mean.  Among  known  problems  are  a  ratner  high  level  of 
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uncertainity  in  the  calibrations,  the  non-stationary  noise  and 
the  possible  error  from  merging  two  seismograms. 

Enhancement  of  transverse  waves 

The  Greenland  Sea  earthquake  of  19  May  1970  has  an  azimuth 
at  QCAZ  of  11.3°  (the  installed  azimuths  of  the  horizontal  strain 
seismographs  are  55°  and  325°)  which  is  such  that  SH  and  Love 
waves  should  be  enhanced  by  the  difference  of  the  orthogonal 
horizontal  strain  seismographs. 

Figure  19  shows  the  sum  and  difference  of  the  horizontal 
strain  seismographs  at  QCAZ  during  the  expected  arrival  time  of 
surface  waves  from  the  Greenland  Sea  earthquake.  The  computed 
arrival  time  of  the  Love  wave  has  been  marked.  Computation  was 
based  on  the  Jef freys-Bullen  travel  times.  The  surface  wave 
arrivals,  as  interpreted,  have  been  enclosed  in  boxes  in  the 
figure. 

Determination  of  local  phase  velocity 

Given  a  single  Rayleigh  mode  whose  horizontal  component  is 
defined  by 


ux  *  g  (knz)  sin  (knx  -  u>t) 


the  horizontal  strain  output  with  an  electromagnetic  transducer 
is 


d  aux 

-  =  wk  g(knz)  sin  ^knx  '  a>t-) 

dt  9x  1 


10- 


and  the  horizontal  pendulum  output  is 


2 

-  (u  )  =  -  a)  g  (k  z )  sin  (k  x  -  wt) 

9  A  II  n 

dt^ 


The  ratio  of  strain  to  pendulum  amplitude  is 

uk  k  1 

n  n 


where  cn  is  the  local  phase  velocity. 

One  can  also  show  using  a  long  algebraic  argument  that  in 
a  homogeneous  half-space,  the  ratio  of  horizontal  strain  amplitude 
to  horizontal  pendulum  amplitude  is  l/cn  for  both  P  and  SV  waves. 

The  application  of  these  formulae  requires  that  only  a 
single  type  of  wave  or  Rayleigh  mode  occur  at  any  given  record 
time.  If  one  had  two  waves  arriving  simultaneously,  such  as  S 
and  PcS  at  A  =  39°,  it  would  not  be  possible  to  determine  phase 
velocity. 

In  computing  local  phase  velocity,  we  identify  the  maximum 
and  minimum  value  of  each  half  cycle  in  the  selected  time 
period,  measure  the  ratio  of  the  strain  and  pendulum  amplitudes 
at  their  maximum  and  minimum  amplitude  values,  and  plot  the 
ratios  against  the  corresponding  half-period  (times  two). 

Figure  20  shows  the  local  phase  velocity  at  QCAZ  calculated 


for  detailed  analysis.  Figure  8  (previously  described)  is  the 
QCAZ-TFO  comparison  for  this  earthquake.  The  analog  computer  was 
used  in  the  present  analysis;  channel  gains  were  adjusted  to 
obtain  cancellation  of  the  sum  or  difference.  Figures  22  to  25 
show  the  outputs  of  various  instruments  for  selected  time 
segments  of  the  Chiapas  signal.  An  analog  circuit  to  obtain 
outputs  phase-shifted  90°  relative  to  each  other  has  been  used 
as  part  of  the  analytical  procedure. 

Figure  22  shows  the  outputs  of  the  radial  long-period  pen¬ 
dulum,  the  vertical  long-period  pendulum  and  their  sum  and 
difference,  for  the  Chiapas  earthquake.  The  calculated  arrival 
times  of  the  phases  P,  S,  and  SS  from  the  Jeffreys-Bullen 
tables  have  been  marked.  There  is  an  additional  phase  at  about 
01  03  *5  that  is  not  predicted  by  the  travel  time  tables.  The 
S  phase  is  90°  out  of  phase  on  the  radial  and  vertical  pendu¬ 
lums.  The  ratio  of  displacement  for  a  homogeneous  half -space 
is  given  by  the  formula 


u. 


u. 


2  s±ri(}>  f y /  (A  +  2p)  -  sin2^]1^ 


z=0 


sin  2<j) 


for  SV  waves  incident  at  an  angle  <|>  (Appendix  II).  If  X  -  y, 

1/2 

the  numerator  becomes  zero  at  35.3°,  where  sin  4>  =  1/ 3 
For  angles  greater  than  <j),  the  ratio  is  imaginary,  indicating 
a  90°  phase  difference  between  the  horizontal  and  vertical 
pendulum  for  incidence  angles  greater  than  35.3°.  For  the 
Chiapas  earthquake,  using  the  slope  of  the  travel  time  curve, 
assuming  a  compressional  velocity  of  6.4  km/sec,  and  X  =  y, 
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one  obtains  <j>  =  34°. 

However  the  many  assumptions  in  this  calculation  render 
the  results  suspect.  For  some  fixed  value  of  $  >  35° 

u  =  j  K  u 

where  j  =  /T  and  K  includes  all  the  constant  terms.  Thus  uz 
(vertical  pendulum)  is  expected  to  lead  u^  (horizontal  pendu¬ 
lum)  by  90°.  This  analysis  of  pendulum  channels  is  included  to 
support  the  data  analysis  of  the  strain-pendulum  combinations. 

Figure  23  shows  the  radial  long-period  pendulum  and  strain 
outputs,  and  their  sum  and  difference.  There  is  an  apparent 
polarity  reversal  in  the  strain  circuit.  For  this  azimuth 
(133.1°)  the  sum  should  enhance  S  and  Rayleigh  waves.  The 
computed  travel  time  of  S  and  the  peak  amplitude  of  the  Rayleigh 
wave  have  been  marked.  There  is  apparent  surface  wave  or  body 
wave  energy  beginning  about  01  11  30,  as  marked  by  an  arrow. 
Prior  to  the  marked  time  the  difference  cancellation  is  not  as 
effective. 

Figure  24  shows  the  vertical  and  horizontal  pendulum  out¬ 
puts  for  the  same  time  segment  as  Figure  23.  An  analog  90° 
phase-shift  circuit  has  been  applied,  so  that  in  the  figure  the 
outputs  will  be  in  phase  or  180°  out  of  phase.  In  this  case  S 
is  enhanced  on  the  difference  and  the  Rayleigh  waves  on  the 
sum  channel.  The  retrograde  motion  beginning  at  01  11  30  is 
nearly  identical  to  that  obtained  from  the  strain-pendulum 
sum  (Figure  23). 

Figure  24  shows  the  transverse  pendulum  output  (P55L) 
and  the  sum  of  the  orthogonal  horizontal  strain  outputs.  The 
latter  is  equivalent  to  a  vertical  strain  instrument.  The 
peak  of  the  Love  waves,  and  an  apparent  SH  component  of  S, 
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have  been  marked  on  the  transverse  pendulum.  As  expected,  the 
Love  wave  at  0111  is  not  seen  on  the  vertical  pendulum  output. 
The  Love  wave  should  also  be  enhanced  on  the  difference  of  the 
orthogonal  strain  instruments,  but  in  this  example  the  signal- 
to-noise  ratio  is  not  actually  improved  over  the  transverse 
pendulum. 

There  is  a  phase  that  is  enhanced  on  the  sum  with  a  peak 
about  01  12.  Enhancement  on  the  sum  trace  indicates  that  it  is 
not  SH-type  motion.  There  is  no  well-defined  evidence  of  its 
occurrence  on  the  vertical  pendulum.  Such  motion  as  is  seen 
there  is  probably  the  early  arriving  Rayleigh  wave. 

The  evidence  is  thus  for  the  occurrence  of  a  low*»frequency 
(T  =  20  secs)  wave  whose  principal  motion  is  horizontal.  It 
cannot  be  clearly  identified  as  a  Rayleigh  or  Love  wave. 


SUMMARY  AND  CONCLUSIONS 


The  beam  comparisons  between  QCAZ  and  TFO  are  consistent 
with  the  hypothesis  that  the  signal  enhancement  of  Rayleigh 
waves  for  the  two  systems  are  approximately  equal  when  system 
noise  is  low.  This  hypothesis  was  investigated  in  detail  for 
one  example  in  which  it  was  possible  to  compare  signal  amplitudes 
to  the  rms  amplitude  of  the  noise. 

Love  wave  enhancement  has  been  obtained  using  a  technique 
that  is  limited  to  signals  whose  azimuth  from  the  station  is 
between  about  30°  to  60°  from  the  installed  azimuths  of  ortho¬ 
gonal  strain  seismographs.  In  the  example  used  both  strain 
channels  had  a  high  noise  level. 

A  Chiapas  earthquake  has  been  examined  in  detail  to 
demonstrate  signal  enhancement  techniques  applied  to  P,  S, 

SH,  Love,  and  Rayleigh  waves. 

Some  preliminary  data  on  use  of  the  strain-pendulum  ratio 
to  determine  local  phase  velocity  has  been  included. 

These  results  show  that  the  strain  seismograph,  in  conjunc¬ 
tion  with  a  pendulum  seismograph,  can  be  used  to  form  beams  at 
a  single  site,  and  that  the  combination  can  be  used  to  enhance 
and  recognize  the  total  seismic  signal. 

At  the  time  this  data  was  obtained,  the  instruments  were 
still  in  the  process  of  being  installed.  Later  data  should 
be  free  of  some  of  the  non-stationary  noise  that  prevented  a 
more  complete  analysis  at  this  time. 
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Figure  1.  Arizona  index  map. 

li 


Figure  2.  Location  of  earthquakes  analyzed 


Figure  3.  Azimuthal  response  of  strain  seismograph  and 
horizontal  pendulum  seismograph  to  longitudinal  plane 
waves  or  Rayleigh  waves. 
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Figure  5.  Response  of  horizontal  strain  seismograph  to 
transverse  waves  (SH  or  Love). 
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Figure  6.  Beam  pattern  comparison  TFO  (6  elements)  to 
QCAZ  (sum  of  strain  and  pendulum) 
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Figure  7.  Beam  pattern  comparison  TFO  (4  elements)  to 
QCAZ  (sum  of  strain  and  pendulum). 


Figure  8.  Example  of  horizontal  rotation  at  TFO. 
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Figure  9.  QCAZ-TFO  comparison  --  near  Coast  of  Chiapas, 
Mexico,  19  May  1970. 
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Figure  11.  QCAZ-TFO  comparison  --  near  Coast  of  Chiapas, 
Mexico,  21  May  1970. 
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FROM  SOUTHEAST  STEERED 
TOWARD  NEW  HEBRIDES  EQ.  (SW) 


Figure  12.  QCAZ-TFO  comparison  --  New  Hebrides  Islands, 
21  May  1970. 


Figure  13.  QCAZ-TFO  comparison  --  Kermadec  Islands 
25  May  1970. 
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Figure  14.'  TFO  back  azimuth  beam,  Kermadec  Islands 
25  May  1970. 
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Figure  15.  QCAZ-TFO  comparison  south  of  Panama,  22  May  1970. 


S55 


Figure  19.  Surface  waves  Greenland  Sea  earthquake,  difference 
of  orthogonal  horizontal  strain  to  enhance  Love  waves. 


Figure  20.  Local  phase  velocity  (Rayleigh)  Kermadec  Islands 
earthquake . 


22.  Analog  computer  analysis  Chiapas  earthquake 


19  May  1970. 


Figure  25«  Analog  computer  analysis  Chiapas  earthquake 
19  May  1970. 


APPENDIX  I 

DATA  FROM  USC§GS  PDE  CARDS 
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Data  From  USC§GS  PDE  cards 
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APPENDIX  II 

HOMOGENEOUS  HALF- SPACE  AMPLITUDE  AND  PHASE  ANALYSIS 


Using  the  formulation  of  Ewing  et  al.  (1957)  it  is  possible 
to  derive  theoretical  formulations  for  the  amplitude  and  phase 
relation  for  the  vertical  and  horizontal  strain  and  pendulum 
seismographs  at  a  free  surface.  This  development  is  limited  to 
P  and  SV  waves . 

In  order  to  make  the  results  somewhat  more  general  the 
assumption  that  X  =  y  has  not  been  made  (X  and  y  are  Lam6 
constants).  The  results  of  Gupta  (1966)  are  a  special  case  of 
results  here  obtained,  since  he  assumed  X  =  y.  Gupta's  analysis 
did  not  include  the  results  for  the  horizontal  strain  and 
pendulum. 

Figures  II- 1  and  1 1 - 2  show  the  geometry  for  incident  P  and 

SV  waves  at  the  free  surface  of  an  elastic  solid.  Using  the 

notation  u  and  u  for  horizontal  and  vertical  displacement 
X  z 

respectively,  and  defining 

2  2  °2 
r  =  tan  e  =*  —  -  1 

2 


one  obtains  the  following  formulation  for  vertical  displacement 
at  the  free  surface  (uz|Q),  vertical  strain  at  the  free  surface 
3u/3z | z_0»  horizontal  displacement  at  the  free  surface  ux|Q,  and 
the  horizontal  strain  3ux/8x|o.  The  algebra  required  becomes 
rather  involved  and  is  not  included.  Additional  notation  is 


Figure  1 1 - 1 .  Reflection  of  P-wave  at  free  surface  of  an 
elastic  solid. 


Figure  1 1 - 2 .  Reflection  of 
elastic  solid. 


SV-wave  at  free  surface  of  an 


identical  to  that  of  Ewing  et  al .  (1957) 


For  P-waves: 


uzl0  =  jkA2 


2r 2  (s2  +  1)  [X (1  +  r2)  +  2m  r2) 


4p  rs  +  (s2  -  1)  [X  (1  +  r2)  +  2\ir2] 


3u 


3z 


k2A, 


4rsX  (1  +  r^) 


z=o 


4yrs  +  (s2  -  1)  [X (1  +  r-2)  +  2yr2] 


(X  +  2m)  (1  ♦  r  ) 


uxl  =  -  jkAj  4rs 


z=o 


4mts  +  (s2-l)  [X  Cl  ♦  r2)  +  2mt2] 


3u 


3x 


(X  ♦  2m)  (1  ♦  r  ) 


=  ■  k  Aj  4rs 


z=o 


4mts  ♦  (s2  -  1)  [X  (1  ♦  r2)  +  2mt2] 


For  SV-waves: 


u  I  =  -  jkB 
z=o  1 


4mts(s  +  1) 


4urs  +  (s2-l)  [X(l  +  r2)  +  2mt2] 


3u 


3  z 


k2Bx  s ( s 2  -  1) 


2X  (1  +  r  ) 


z=o 


4mts  +  (s2  -  1 )  [X  (1  +  r2)  +  2mt2] 


ih 


1 


u 


s(s2-l)  (X  +  2y)  (1  +  r2) 

=  -  2kj  Bj - 

2=0  4prs  +  (s2  -  1)  [X (I  +  r2)  +  2yr2] 


3u_ 


3x 


=  -  2k2B1  s (s 2  -  1) 


(X  +  2y )  (1  +  p4) 


z=o 


4yrs  +  (s2  -  1)  [(1  +  r2)  +  2yr2] 


The  time  varying  factor  exp  [jk(ct  -  x)]  ha?  been  dropped  in 
each  case.  Using  the  above  formulae  and  letting  0  be  the  incident 
angle  for  P  waves  (Figure  1 1 - 1 )  one  obtains,  after  some  algebraic 
manipulation: 


u  I  COS0 


U”  '  z=o  sin20 


(X  +  2y  cos40) 

„  rX  +  2W  .  2q11/2 
2y[ -  -  sm  0] 


3uz 

Tz~ 

u_ 


^  +  2y  2  I/? 

2jkXysin0  [ -  -  sin  0]A 


(X  +  2y)  (X  +  2y  cos20 ) 


z=o 


X 


3u. 

I 

dz 

u„ 


z  jk 


z-o 


X  +  2y 


9u 


dx 


u 


=  jk 


X  +  2y 


z=o 


Similarly  for  SV  waves,  where  0  is  angle  of  incidence  (Figure  II-2) 
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